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Abstract 


Necessary  and  Sufficient  Conditions  for  the  Existence  of  Complementary 
Solutions  and  Characterizations  of  the  Matrix  Classes  Q  and 


Faiz  A.  Al-Khayyal 
Georgia  Institute  of  Technology 


The  linear  complementarity  problem  has  a  solution  if  and  only  if 
there  exists  a  right-hand-side  vector  of  an  associated  dual  linear  pro¬ 
gram  such  that  at  least  one  LP  solution  falls  inside  a  specific  open  set 
This  condition  is  used  to  give  different  characterizatons  of  the  matrix 
classes  Q  and  Q^.  Moreover,  they  are  used  to  derive  verifiable  suffi¬ 
cient  conditions  for  a  subclass  of  Q. 
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Necessary  and  Sufficient  Conditions  for  the  Existence  of  Complementary 
Solutions  and  Characterizations  of  the  Matrix  Classes  Q  and 

Faiz  A.  Al-Khayyal 
Georgia  Institute  of  Technology 


1.  Introduction 


The  linear  complementarity  problem  is  to  find  a  real  n-vector  x  such 


that 


Mx+q  >  0,  X  >  0,  x^(Mx+q)  =  0,  (1) 

where  q  Is  a  known  real  n-vector  and  M  is  a  known  real  nxn  matrix.  The 
problem  Is  generally  denoted  by  the  pair  (q,M).  For  given  M,  the  feasi¬ 
ble  region  of  (q,M)  Is  given  by  the  convex  polyhedron 

S(q,M)  »  {xeR^:  Mx+q  >  0,  x  >  0).  (2) 

When  S(q,M)  *  0,  the  problem  (q,M)  and  the  vector  q  are  said  to  be  feasi¬ 
ble  •  A  feasible  problem  (q,M)  does  not  necessarily  have  a  solution.  The 
main  theorem  of  this  paper  provides  a  necessary  and  sufficient  condition 
for  problem  (q,M)  to  have  a  solution.  This  result  is  then  used  to 
characterize  two  Important  matrix  classes. 

Numerous  procedures  have  been  proposed  for  solving  problem  (q,M), 
Including  Lemke  [11],  Cottle  and  Dantzig  [3],  Murty  [14],  Eaves  [7], 
Saigal  [18],  Garcia  [8],  and  Karamardian  [9].  More  recently,  Mangasarlan 
[12,13],  Cottle  and  Pang  [4],  Shlau  [19],  Pardalos  and  Rosen  [16], 

Ramarao  and  Shetty  [17],  and  Al-Khayyal  [1]  have  developed  different 
approaches  for  solving  either  (1)  or  an  equivalent  restatement  of  the 


problem.  While  the  preceding  list  Is  not  exhaustive,  with  the  exception 
of  [1,16,17],  all  the  methods  assume  structure  In  the  matrix  M.  This 
note  addresses  two  such  classes  of  structured  matrices. 

The  classes  of  matrices  considered  are 

Q  =  qeR*'  (q»M)  has  a  solution} 

and 

Qq  “  {MeR"*":  qeR"  and  S(q,M)  #  0  =*>  (q,M)  has  a  solution}. 

The  class  has  also  been  denoted  K  In  the  literature  (e.g..  Pang  [15]). 
Several  authors  have  Investigated  these  classes  and  have  obtained  full  or 
partial  characterizations  of  such  matrices  (see,  e.g.,  Mangasarlan  [12], 
Cottle  and  Pang  [4],  Pang  [13],  Stone  [20],  Shiau  [19],  Doversplke  and 
Lemke  [6],  Kelly  and  Watson  [10],  and  Cottle,  Randow  and  Stone  [5]). 

A  geometric  characterization  of  the  members  of  Is  that  the  union 

of  complementary  cones  of  any  matrix  In  Is  a  convex  set.  The  class  Q 
Is  a  subset  of  for  which  the  union  of  complementary  cones  of  any 
member  In  Q  Is  r” .  In  this  note,  we  obtain  new  necessary  and  sufficient 
conditions  for  the  classes  Q  and  Qq  which  have  a  different  geometric 
Interpretation.  Verifiable  conditions  based  on  linear  programming  are 
provided  which  are  sufficient  for  membership  In  Q. 

2.  The  Conditions 

Let  Y  denote  the  unit  hypercube 


where  e  =  (1,1,..., 1)  is  the  summation  vector.  For  fixed  yeY,  define 
the  polyhedral  set 

V(y,M)  =  {veR*':  M^v  <  e-y,  v  >  -y}  , 
and  consider  the  linear  program 

min{q^v;  veV(y,M)}.  (3) 

Denote  the  solution  set  of  this  problem  by 

T 

argmin{q  v:  veV(y,M)}  , 

and  consider  the  following  open  polyhedral  set 

V  -  {vcR”:  (M-I)^v<e}.  (4) 

Theorem  1.  Let  qeR*'  and  Mcr"*".  The  problem  (q,M)  has  a  solution  if  and 
only  if  there  exists  a  vector  yeY  such  that 

argrain{q  v:  veV(y,M)}  n  V  *  0.  (5) 

Proof ;  (Necessity).  Let  x  solve  the  problem  (q,M).  Let  yeY  be 


defined  as  follows 


E 

S 

i 


is 


I 


I ». 


V 
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if  Xj^>0 
if  Xj^»0. 


It  follows  from  linear  programming  duality  that 

T  -  “ 

0  e  argmin{q  v:  veV(y,M)}  n  v 


so  that  (S)  is  satisfied. 

(Sufficiency).  Let  yeY  be  such  that  (5)  holds  and  let 
-  T  -  * 

V  e  argmln{q  v:  vcV(y,M)}  V.  Then  by  the  Karush-Kuhn- 
Tucker  Conditions  there  exist  multipliers  (X,£)  with  XeS(q,M) 
and  5  =*  MX+q  such  that 


X^le-y-M^v)  -  0 


(6) 


T  T 

It  follows  that  X  (MX+q)  =  X  ^  =  0  so  that  X  solves  (q,M). 
For  if  this  were  not  true,  then  there  would  exist  an  index  i 
such  that  =  Xj^(MX+q)j^  >  0.  By  the  complementary  slack¬ 

ness  condition  (6),  we  must  have 


0  and  [y  +  v)j^  »  0. 


Adding  gives  [e  +  (l-M)^v) 


0  which  contradicts  the  assump- 


tion  that  veV. 


Note  that,  when  (5)  holds  for  some  ycY,  problem  (q,M)  has  a  solution 
and  q  Is  feasible.  Thus,  any  Infeasible  q  will  violate  (5)  for  all  yeY. 
The  theorem  can  be  restated  to  characterize  the  classes  Q  and  Q^. 
Corollary  1.  Let  qeR^  and  MeR.”*”.  The  matrix  MeQ  if  and  only  if  for 
each  q  there  exists  a  vector  yeY  such  that  (5)  holds. 


Corollary  2.  Let  qeR*^  and  The  matrix  ^eQq  if  and  only  if  for 

each  feasible  q  there  exists  a  vector  yeY  such  that  (5) 
holds . 


From  the  proof  of  the  theorem,  we  see  that  when  a  complementary 
solution  exists,  there  exists  a  yeY  such  that  the  origin  must  necessarily 
solve  the  linear  program  (3).  The  following  example  illustrates  that 
this  need  not  be  the  case  when  q  is  feasible  and  no  complementary 
solution  exists. 


Example  1 


1/2  1  I 

^  M  =  ,  q  = 

.  1  -1/2J  L-1. 


Here,  we  have  S(q,M)  =  {x:  l/Zx^^+x^+l  >  0,  Xj^-l/2x2-l  >  0, 

Xj^  >  0,  X2  >  0}  #  0,  but  there  is  no  yeY  for  which  the  origin 
solves  problem  (3).  In  fact,  the  optimal  solution  to  the 
linear  program  (3)  Is  not  in  V  for  any  yeY. 


Formalizing  these  observations,  let  cl  V  denote  the  closure  of  the 

A  A  ^ 

set  V  and  let  3(cl  V)  denote  the  boundary  of  the  set  cl  V.  We  have  the 


•  .-.‘.I 

-•.-.‘.••■.I 

m  I 
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following  consequence  of  the  theorem. 


Corollary  3.  Let  qeR*^  and  MeR^  The  problem  (q,M)  does  not  have  a 
solution  If  and  only  If  for  every  y  e  Y, 

argmin{q^v:  veV(y,M)}  c  3(cl  V).  (7) 

Proof;  From  Theorem  1,  It  follows  that  problem  (q,M)  does  not  have 
a  solution  if  and  only  if  for  every  yeY 

T  -  “ 

argmin{q  v:  veV(y,M)}  n  V  =  0.  (8) 

However,  it  is  clear  that 

T  - 

argmin{q  v:  veV(y,M)}  ^  cl  V.  (9) 

The  result  follows  since  (8)  and  (9)  are  true  if  and  only 
if  (7)  holds. 

For  completeness,  this  will  be  stated  in  terms  of  the  classes  Q  and 

Corollary  A.  Let  qeR*^  and  MeR*'’'''.  The  matrix  M^Q  if  and  only  if  there 
exists  a  q  such  that  (7)  holds  for  every  yeY. 

Corollary  5.  Let  qeR*'  and  MeR*^  The  matrix  if  and  only  if  there 


exists  a  feasible  q  such  that  (7)  holds  for  every  ycY. 

New  sufficient  conditions  for  membership  in  Q  can  also  be  derived 
from  Theorem  1.  Let 

V'  -  {vcR":  M^v  <  e,  V  >  -e}.  (10) 

Corollary  6.  If  V'  c  V  then  MeQ» 

Proof :  For  every  yeY,  we  have  V(y,M)  ^  V'.  Hence,  for  any  qeP”, 

“  T 

if  V'  c  V  then  argmln{q  v:  veV(y,M)}  ^  V  so  that  (5)  is 
satisfied.  Thus  every  qeR^  yields  a  solution  to  problem 
(q,M)  so  that  MeQ. 

The  condition  in  Corollary  6  is  rather  strong  since  it  Implies  that 
(5)  is  satisfied  for  all  yeY.  Nevertheless,  it  Includes  some  P-matrices 
as  Illustrated  by  the  following  example. 

[1  1/2' 

, 

-1/2  1 

V'  =  {(v^,V2):  ^Vj^-v^  <  2,  v^+2v2  <  2,  Vj^  >  -1,  v^  >  -1} 

*=  {(Vj^,V2):  Vj^  <  2,  v^  >  -2}  =  V. 

On  the  other  hand,  the  next  example  shows  that  there  are  P-matrices 
that  violate  the  condition  V'  c  v. 


’1/2  1/2 

Example  3.  M  = 

[-2  1 


V'  =  {(Vj^,V2):  Vj^-4v2  <  2,  Vj^+2v2  <  2,  >  -1,  >  -1}, 

V  =  {(V2,V2):  +  4v2  >  -2,  <  2}. 

Here,  V'  /  V  since  (-l,-3/4)  is  in  V'  but  not  in  V. 

The  preceding  corollary  suggests  checkable  conditions  to  test 
whether  or  not  V'  V.  First,  recall  that  a  nonzero  vector  d  is  called  a 
direction  of  a  set  if  the  ray  emanating  from  every  point  in  the  set  in 
the  direction  d  is  entirely  contained  in  the  set  (see,  e.g.,  Bazaraa  and 
Jarvis  [2],  Chapter  2).  Let  D(V')  and  D(V)  denote  the  set  of  directions 

A 

of  V'  and  V,  respectively.  These  sets  define  the  polyhedral  cones 

D(V')  =  {d  *  0:  M^d  <  0,  d  >  0} 

and 

D(V)  =  {d  It  0:  (M-I)^d  <  0}. 


Note  that 


D(V')  c  D(V). 


(11) 


Theorem  2.  The  set  V'  is  a  subset  of  V  If  and  only  if  for  every  extreme 


point  V  of  V 


(12 


Proof : 


either  (M  v)^  =  1  and  >  0 

or  (M^v)^  <  0  and  =  -1 

for  all  1  ~  l}2^*«*^n« 

(Necessity).  Suppose  V' ^  V  and  let  v  denote  an  arbitrary 
extreme  point  of  V'.  Since  v  is  an  extreme  point,  then  at 
least  n  linearly  independent  constraints  of  the  2n 
constraints  of  V'  must  be  binding  at  v;  thus, 

either  (M^v)^  =1  or  =  -1  (13' 

A 

for  all  i  =  1,2, ...,n.  But  V'  c  V  so  that 

^^^^^1  -  1  (1^ 

for  all  1  =  l,2,««.,n.  Clearly,  (13)  and  (14)  imply  (12). 

(Sufficiency).  Assume  that  the  set  V'  has  k  extreme  points 

12k  1  2  P 

{v  ,v  ,...,v  }  and  I  extreme  directions  {d^  ,d^ , . . . ,d^} .  By 

(12),  we  have 


(M-I)^v^  <  e 


(15 


for  all  j  -  l,2,...,k.  Since  V'  is  a  convex  polyhedral  set 

then  veV'  if  and  only  if  there  exist  nonnegative  vectors 
T 

(X,vi)  with  e  X  =  1  such  that 


fj 


f 


r 


I 


k  • 

I V' 


t: 


£ 


y. 

•V 
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k  .  £  . 

V  =  £  X  .v^  +  £  u  .d^  . 

j=l  ^  j=l  ^ 


Hence , 

Ic  ■ 

(M-I)^v  =  £  X.(M-I)^v^  +  £  u.(M-I)^d^  <  e. 

j=l  ^  j“l  ^ 

The  last  Inequality  follows  from  (15)  and  (11)  and  the  condi- 

T  * 

tions  e  X  =  1  and  u  >  0*  Therefore  veV,  which  completes  the 

proof • 

3.  Summary  and  Conclusions 

This  note  presents  new  necessary  and  sufficient  conditions  for  the 
existence  of  a  solution  to  the  linear  complementarity  problem.  These 
translate  to  similar  conditions  for  the  classes  of  matrices  Q  and  and 
for  their  complements.  In  particular.  Theorem  1  (Corollary  2)  states 
that  existence  of  a  solution  to  the  linear  complementarity  problem  (mem¬ 
bership  in  Qq)  Is  equivalent  to  the  existence  of  a  solution  to  the  linear 
program  (3),  for  some  fixed  yeY,  which  Is  also  In  the  open  set  V  given  by 
(4).  From  this  follows  the  interesting  condition  on  membership  in  the 
complement  of  Q^,  which  Is  equivalent  to  the  existence  of  a  feasible  q 
such  that  every  solution  of  the  linear  program  (3),  for  every  fixed  yeY, 

A 

Is  on  the  boundary  of  the  closure  of  V. 

While  these  conditions  are  not  easy  to  verify,  they  do  provide  a  new 
framework  from  which  new  results  can  be  derived.  For  example.  Corollary 
6  provides  a  sufficient  condition  for  membership  in  Q.  This  condition 
can  be  verified  by  using  the  test  suggested  by  Theorem  2.  This  test 
Involves  generating  all  the  extreme  points  of  V'  and  checking  each 


'  s  s'  N 
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against  (12).  Note  that  to  establish  that  a  matrix  is  not  in  Q  would 
require  generating  all  the  extreme  points  only  in  the  worst  case.  How¬ 
ever,  to  confirm  that  a  matrix  is  in  Q,  all  extreme  points  must  be  gene¬ 
rated. 
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